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ABSTRACT 

Using a wavelet decomposition technique, we have extracted the Hurst exponent for a 
sample of 46 long and 22 short Gamma-ray bursts (GRBs) detected by the Gamma- 
ray Burst Monitor (GBM) aboard the Fermi satellite. This exponent is a scaling 
parameter that provides a measure of long-range behavior in a time series. The mean 
Hurst exponent for the short GRBs is significantly smaller than that for the long 
GRBs. The separation may serve as an unbiased criterion for distinguishing short and 
long GRBs. 
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1 INTRODUCTION 

Our present understanding of complex astrophysical objects such as cataclysmic variables (CVs), active galactic nuclei (AGNs), 
and gamma-ray bursts (GRBs) comes nearly entirely from the temporal and spectral analyses of their photoemissions (with 
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some additional information coming from possible associations, such as host galaxies). In this paper we consider the temporal 
aspects of GRB light curves observed by the Gamma-Ray Burst Moni tor (GBM) aboard the Fermi satellit e . Many studies of the 



tempo r al properties of GRB ligh t curves have bee n published, such as 



Nemirofl 



(2009): 



Hakkila fc Preecel |201ll ) 



Nemirofl 



(2000): 



Norris et al 



(2005): 



H akkila fc Nemirofl 



(|2012f ) , from the perspective that light curves are comprised of a series of displaced 
pulses and that by fitting the individual pulses and associating pulses at various photon energies one can arrive at a holistic 
understanding of light curves which in turn may be used to constrain the physics of the engines that produce them. The 
main appeal of this approach is the intuitive connection between pulses and collisions in the internal shock model. While this 
is a perfectly reasonable method, issues do arise concerning the functional form to use for pulse fitting and how to discern 
actual pulses from stochastic fluctuations in the light curves. The situation is further exacerbated by the fact that GRB light 
curves exhibit considerable variation in duratio n and in pulse profil e. We note the significant progress made in non-parametric 



analyses using the Bayesian block technique (jScargle et al 



20121 ). In this type of analyis the duration of a light curve is 



represented as a tessellated block of time which can be partitioned into a complete array of sub-blocks in any number of ways. 
An optimal partition (a partition of sub-blocks maximizing a fitness function) is shown to exist, be unique, and computable 
iteratively. The optimal partition of sub-blocks is determined, given a prior probablility distribution for the number of blocks, 
by finding the model best representing the data as sets of piece-wise constant segments or sub-blocks. This technique shows 
great promise in resolving statistically significant temporal features from noise and detector related artifacts. 

An ideal complementary approach to probing light curves would be one which handles seemingly disparate profiles on 
an equal footing and distills their complex forms into a single parameter which may be used to compare one light curve with 
another. One such method was pioneered by Harold Edwin Hurst in 1951 with a technique he invented called the rescaled 
range analysis (R/S) which was later improved upon by Benoit Mandelbrot. The eponymous parameter resulting from the 
rescaled range analysis is called the Hurst exponent, H, and is closely related to the fractal dimension, D, the understanding of 
which Mandelbrot spent much of his career developing. In fact, fractional Brownian motion (fBm), which Mandelbrot defined 
in 1968, is parametrized solely by H and serves as a useful model for discussing time series. After determining H for a given 
time series one is in a position to make several statements about the nature of that time series including whether the sequence 
appears random or whether it is persistent or anti-persistent, and if it is, whether it exhibits long-range dependence, and 
over what time scales these characteristics are operative. All of these are informative quantitative statements, especially if the 
specific process generating the time series is partially or completely unknown, in which case, these statements are perhaps all 
one can really say about the process given the available information. Some fields of research in which interesting work is being 
done with Hurst exponents are financial markets, seismology, anesthesiology, astrophysics, plasma physics and genomics. 
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Pseudo Time 

Figure 1. Simulated fractional Brownian motions with different values of H: a) H = 0.25, b) H = 0.50, c) H = 0.75. 



We point out that neither the pulse fitting methods nor the Bayesian block analysis (|Scargle et al 



20121 ) yields information 



directly relatable to the Hurst exponent as does the wavelet analysis. One approach to access the Hurs t exponent from a 



Bayesian block framework that seems reasonable would be an adaptation of the Box-Counting algorithm (jFeder 
a Bayesian-Box-Counting algorithm is outside the scope of this paper. 



19881 ). Such 



The estimation of the Hurst exponent and the related scaling expone nt, a, has a history in astrophysics ijAnzolin et al 



2010 



Tamburini et al 



2009 



Walker fc Schaefer 



2000 



Fritz fc Brucb 



19981 ) for both Cataclysmic Variables (CVs) and GRBs. 



We propose that a similar determination of H for GRB light curves will be a valuable tool for categorization and we present 
a separation of long and short GRBs based on H . 



2 METHODOLOGY 



2.1 Hurst Exponent and Self- Affinity 



Pioneering work i n self-simila rity and long-range dependence was first published in 1951 by Hurst in the study of annual 



Nile River levels, ( Hurst 



19511 ). Hurst examined several decades of data to determine what should be the minimum size of a 
reservoir so that it neither overflows nor runs dry due to yearly fluctuations and made the unexpected observation that annual 
Nile River levels were not independent from one another but instead exhibited a memory of past events. 

In this analysis of time-series data we search for statistical fractals, i.e., fractals whose statistical characteristics are 
independent of time scale. Such fractal time-series are called self-similar. There is another class of statistical fractals whose 
scale invariance is broken but can be restored by a multiplicative factor. These statistical fractals are called self-affine. 



Mandelbrot! (| 19851 ) defined a time-series, X(i) with t £ {to . . . tN-i}, to be self-affine if, after a rescaling t — > Xt the following 



4 MacLachlan et al. 

relation is satisfied, 



X(t) = \~ X(\t). 



(1) 



The exponent, H, is the Hurst exponent, (|Hurs1 



example of a self-affine time-series, also given by 



1951 ) and t he sy mbol = denotes equality in distribution. The canonical 



Mandelbrot! (|l968T ). is fractional Brownian motion, fBm. Stationary in the 



context of this paper is second-order stationarity which means the first and second moments obey the following relations 
E{X(t)} = fi x 

E{X(t 2 )X(tt)} = 7«a - *i) = 7(r), (2) 

where fix is the sample-mean, 7 is the auto-covariance sequence and r = £2— £1 is the lag. The Hurst exponent, H, parametrizes 
the degree of statistical self-similarity which a time-series exhibits. A self-similar series may be sub-divided into three categories: 
A series with 1/2 < H < 1 is r eferred to as persistent or long-range dependent while a series with < H < 1/2 is referred to 
as anti-persistent, 



Federl (|1988f l. For H — 1/2 we have neither persistence nor anti-persistence and this corresponds to the case 
of random and uncorrelated events. The Hurst exponent provides a model-independent characterization of the data. Three 
examples of times series with different values of H are shown in Fig. [T] A graphical depiction of the rescaling described by 



Eq. [T]for a time series with H — 0.25 is given in Fig. [2] and for a time series with H = 0.75 in Fig. [3] 



2.2 Wavelet Transforms 



Wavelet transformations hav e been shown to be a natural tool for multi-resolution analysis of non-stationary time-series 



( Flandrin 



1992 



Mallat 



1989). Wavelet analysis is similar to Fourier analysis in many respects but differs in that a wavelet 
basis function, ip(t), is well- localized while Fourier basis functions are global. Locali zation means th at outside some range the 
amplitudes of wavelet basis functions go to zero or are otherwise negligibly small, 



Percival 



(2000). On the other hand, the 



wavelet transform is similar to the Fourier transform because they both are expansions into a complete orthogonal basis and 
resolve low-frequency, large scale structure from high-frequency, small scale structure. 

Wavelet analysis is said to be multiresolution because the time-series under investigation is interrogated at multiple 
scales by a basis set of wavelets which are rescaled and translated versions of an original wavelet commonly referred to as the 
mother- wavelet, ip(t), 



Ja \ a 



(3) 



where a represents an octave or time-scale and the parameter b gives the position of the wavelet within the octave. 
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The continuous wavelet transform (CWT) coefficient, C a ,b, of a time-series for some scale and position is computed as 



Ca.b = / X{t)ll) a , b {t)dt. 

la 



(4) 



2.2.0.1 Wavelet Analysis The wavelet-transform technique for estimating self-affinity is outlined here. By substituting 
the distribution relation in Eq. [T]into Eq. [4] we find 



4= / X(t)i/j a , b (t)dt 
Vo J 



x -(H+l/2) 

\J Xa 



X{Xt)^ a , b ( 



Xt-Xb 
Xa 



)d(Xt) 



(5) 



It is straightforward to see from Eq.[S]that a self-affine time-series will have wavelet coefficients whose variance over a particular 
scale, Xa, is related to the scale parameter A by, 



log va,r(C\a t \ b ) = (2H + 1) log A + constant. 



(6) 



2.2.0.2 Fast Wavelet Transforms Similar to the CWT, the discrete fast wavelet transform (FWT) is also a multiresolu- 
tion operation owing to the construction of the wavelets, tpj,k, which form the basis of the discrete fast wavelet transform. We 
employed the discrete wavelet transform because of its high degree of computational efficiency. In order to distinguish between 
the CWT and its FWT counterpart we make a slight change of notation. Just as before, the Vj,k> are rescaled, translated 
versions of the mother wavelet, ip, 



& > = 2- j7 V(2" , 't-*). 



The coefficients of the FWT are written as 



dj, fe = {X,ij)j,k}, 



where j and fc play the roles of a and ft, respectivel y. Moreover, the values which j and fc assume obey the dyadic partitioning 



scheme (Mallat 



1989 



Addison 



2002 



Perciya] 



20001 ) . That is, for a time series whose number of elements is given by iV = 2" 



j < m - 1, 



and 



< k < 2° - 1. 
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Applying the dyadic partitioning scheme removes any redundant encoding of information by the wavelet transform coefficients 
and guarantees orthogonality among the wavelet basis for any change in j or k, 

(i>j,k,ipj',k') = 8jj'5 k ,k'- (8) 



2.3 Logscale Diagrams 



The average power of the light curve at time scale j is expressed as /3j and may be written in terms of the variance of the 
FWT coefficients as 



n~ — l 
1 - 

P 3 = var(Uj, fc ) = — 2^ 



(9) 



where rij is the number of coefficients at scale, j, (|Abrv et al 



2000 



2003T ). Similarly to Eq. [6] it has been shown by 



Flandrin 



1992J) that for a series with non-stationary statistics the power-law variance of wavelet coefficients goes like 



l°g2(A) = (2-ff + l)j + constant, 



(10) 



where H is the Hurst exponent. 



Masrvl |1993|) later extended this result to a larger class of non-stationary problems with 
stationary increments in the low-frequency limit and showed that fBms are a special case. A plot of Eq. [10] is referred to as a 
logscale diagram. Logscale diagrams are useful for identifying scaling regions, i.e., the range of octaves over which self-affine 



scaling occurs. The slope, a, of the scaling region determines related to the Hurst exponent through a = 2H + 1. 



2.4 Choice of Wavelet Basis 



As in any orthogonal transformation, the basis functions to use in a wavelet transform is a matter of strategic choice. One 
typically chooses a basis that e mphasizes som e characteristic of interest. Commonly used families of wavelet bases are the 



20021 ). We chose the Haar wavelet basis which is the simplest of the Daubechies 



Coiflet, Daubechies, and Haar (jAddison 
family. 

The Haar wavelet basis was chose n from among all other possible bases because it has the fewest number of vanishing 



moments and most compact 



variance (IXizheng fc 2 



block method (jScargk 



1998; 



suppo rt (jAddison 



2002J), has a straightforward interpretation, i.e., is equivalent to the Allan 



1997) and is const ant over its interval of support similar to the model assumed in the Bayesian 



Scargle et al 



2013 1. 



The Haar basis is not without some defects, as noted by 



Kaplan fc Jav Kuo 



199J) and 



Flandrin! |l992j). Namely, the 



Haar wavelet transformation is known to underestimate the actual Hurst exponent and this phenomenon is a function of the 
coarseness of the binning, the number of counts in the light curve, and also of H itself. We show in Sec. 12.61 that this effect 
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is present but smaller than «1<t for a set of simulated light curves and is likely to be smaller for actual data. However, we 
consider that the advantages of the Haar basis outweigh its disadvantages. 



2.5 Minimizing Uncertainties 



2. 5. 1 Circular Permutation 



2000 



Coifman 



19951 ) of light 



Spurious artifacts due to incidental symmetries resulting from accidental misalignment (|Percivall ll 
curves with wavelet basis functions are minimized by circularly shif ting the light curve against the basis functions. Circular 



shifting is a form of translation invariant de-noising (|Coifman 



19951 ). It is possible a shift will introduce additional artifacts 



by moving a different symmetry into a susceptible location. The best approach is to circulate the signal through all possible 
values, or at least a representative sampling, and then take an average over the cases which minimizes the effect of spurious 
correlations. 



2.5.2 Reverse-Tail Concatenation 

Both discrete Fourier and discrete wavelet transformations imply that the expansion is periodic, with the longest period equal 
to the full time range of the input data. This can be interpreted to mean that for a series of N elements, {Xo, X\ . . . Xjv-i} 
then Xo is made a surrogate for Xm and Xi is made a surrogate for Xjv+i, and so forth. This assumption may lead to trouble 
if Xo is much different from Xjv-i. In this case, artificially large variances may be computed. Reverse-tail concatenation 
minimizes this problem by making a copy of the series which is then reversed and concatenated onto the end of the original 
series resulting in a new series with a length twice that of the original. Instead of matching boundary conditions like, 

Xo, Xi, . . . , Xn-i, Xo, (11) 

we match boundaries as, 

Xo, Xi, . . . Xjv_i, Xjv-i, . . . , Xi, Xo- (12) 

Note that the series length has thus artificially been increased to 2X by reversing and doubling of the original series. Conse- 
quently, the wavelet variances at the largest scale in a logscale diagram reflect this redundancy. This is the reason that the 
wavelet variances at the largest scale are excluded from least-squares fits of the scaling region. 
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H=0.25 




Pseudo Time 



Figure 2. Graph of Bjj{t) with H = 0.25. A box is placed around a sub-range of t (lower left hand corner). The box is zoomed into 
with time axis scaled by a and amplitude scaled by a~ H . This is a self-affine transformation that not only makes the rescaled version 
qualitatively 'similar' to the original but also preserves the variance as computed in Eq. [9] 

2.5.3 Poisson Operator 

Photon counting statistics are considered in a bootstrapping procedure by applying a Poisson operator, V(Xi, Xi), to every 
light curve prior to analyzing. Each light curve is binned initially at 200 ^-seconds and the number of counts per bin, Xi, 
is used as a mean value, A;, to be supplied to a Poisson random number generator. The value returned from V(\i,Xi) is 
used to replace the number of counts stored in Xi. The Poisson operator is applied to the signal Xi prior to every circular 
permutation. We show in Sec. 12.61 that the Poisson operator does not affect the measured slope of logscale diagrams above 
the Poisson level. 



2.6 A Test Case: Fractional Brownian Motion 



Spatial-temporal fractional Browni an motions (fBm's) are a useful model for studying self-similarity and long-range depen- 
dence in non-stationary time-series, 



Mandelbrot! (|1968f ) and are characterized by a single parameter, H, the Hurst exponent. 
An fBm with a particular H is expressed as Bn{t) and has the property of self-similarity over a range of scales after a rescaling 
of axes, 

B H (t) = a- H B H (at), (13) 

where = denotes distributional equality as in Section 12.11 The efficacy of the H estimation procedure was tested using 
simulated data in the form of fractional Brownian motion (fBm) time series. Two tests were performed; in the first test we 
examine the ability of our algorithm to determine H from fBms in the presence of Poisson noise and in the second test we 
examine how well we can determine H at H = 0.25, H — 0.50, and H — 0.75 from noise-free fBms. 
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H=0.75 




Pseudo Time 



Figure 3. Graph of Bjj{t) with H = 0.75. A box is placed around a sub-range of t (lower left hand corner). The box is zoomed into 
with time axis scaled by a and amplitude scaled by a~ H . This is a self-affine transformation that not only makes the rescaled version 
qualitatively 'similar' to the original but also preserves the variance as computed in Eq. [9] 




Pseudo Time j [octave] 

Figure 4. Panel a) shows a sample fBm pre-processed and ready to be analyzed in black and the same light curve after applying the 
Poisson operator, V, in red. Panel b) shows the Poisson noise that has been added by V . In panel c) logscalc diagrams illustrate the effect 
of Poisson statistics on the Hurst exponent. The bare fBm is shown in black, the dressed Poisson-type fBm is in red, and the residual 
Poisson noise in shown in blue. 



2.6.1 Test 1 



The numerical computing environment MATLAB was used to produce 1000 realizations of fBms with scaling parameter H 
randomly chosen from the range 0.0 < a < 1.0 by using a uniform random number generator. Copies of the fBms were 
combined with a Poisson operator as described in Sec. 12.5.31 The fBms and the Poissonian fBms thus produced are shown in 
black and red respectively in panel a) of Fig. [4] Panel b) shows the Poisson noise that has been added by V. The Logscale 
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Figure 5. Histograms of 3000 simulated fBm traces. Three categories of fBms were generated with known Hurst exponents, H = 0.25, 
H = 0.50, and H = 0.75. These fBms were then analyzed to recover the Hurst exponent. The histograms are the results of this analysis. 
Pairs of vertical lines are drawn for each peak. The shorter of the two indicates the known H used to generate the fBms and the longer 
of the two indicates the H extracted by our analysis. Results are tabulated in Table [T] 



diagrams in panel c) illustrate the effect of Poisson statistics on the Hurst exponent. The bare fBm is shown in black, the 
dressed Poisson-type fBm is in red, and the residual Poisson noise in shown in blue. The logscale diagram for the bare fBm in 
panel c) exhibits a clean slope across all octaves. We see the effect of a Poisson noise operator; it adds to the signal variance, 
constant across all octaves. Below some octave the signal is completely dominated by noise but above that octave the slope 
of the logscale diagrams is independent of V ■ See for example the black and red symbols for j ^ 6. 



2.6.2 Test 2 

In the second test, 3000 simulated Poisson-type light curves were generated. The simulated data were divided into three 
subgroups of 1000 according to H. The three subgroups were H — {0.25, 0.50, 0.75}. The simulated data in each group were 
analysed and an attempt was made to recover the value of the Hurst exponent, H, used to generate the fBm. The Hurst 
exponent was estimated by a least squares fit to the scaling portion of the logscale diagrams to determine a and then H is 
found from Eq. [TO] 

Results of the second test can be seen in Fig. [5] and Table [T] The results show that the FWT analysis with the Haar 
wavelet basis does underestimate the value of H as discussed in Sec. 12.41 but the magnitude of the error is not significant for 
our purpose. 



3 DATA REDUCTION 

The Gamma-Ray Burst Monitor (GBM) on board Fermi observes GRBs in the energy range 8 keV to 40 MeV. The GBM 
is composed of 12 thallium-activated sodium iodide (Nal) scintillation detectors (12.7 cm in diameter by 1.27 cm thick) that 
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Table 1. Summary of results in Fig. [5] 
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Figure 6. Histogram of H extracted from long and short GRBs. The result for long GRBs is plotted as the solid line while the short 
GRB result is plotted with the dashed line. Note the overlap but also that the means are displaced from one another. 



are sensitive to energies in the range of 8 keV to 1 MeV, and two bismuth germanate (BGO) scintillation detectors (12.7 cm 
diameter by 12.7 cm thick) with energy coverage be tween 200 keV and 40 MeV. The GBM detectors are arranged in such a 



way that they provide a significant view of the sky (jMeegan et al 



2009) 



In this work, we have extracted light curves for the GBM Nal detectors over the entire energy range (8 keV - 1 MeV, also 
including the overflow beyond 1 MeV). Typically, the brightest three Nal detectors were chosen for the extraction. Lightcurves 
for both long and short GRBs were extracted at a time binning of 200 microseconds. The long GRBs were extracted over 
a duration starting from 20 seconds before the trigger and up to about 50 seconds after the Tgn for the burst without any 
background subtraction. For short GRBs, durations were chosen to be 20 se conds before the trig ger and 10 seconds after the 



Tqq. The Tqo durations were obtained from the Fermi GBM-Burst Catalog (jPaciesas et al 
and 22 short GRBs used in this study are tabulated in Tables [2] and [3] 



2012J). Summaries for the 46 long 



4 RESULTS AND DISCUSSION 

We have used a technique based on wavelets to determine the Hurst exponents for a sample of GRB prompt-emission light 
curves. As noted in Section 12.11 the Hurst exponent provides a measure of correlated (or long-range) behavior in a time 
series. The extreme values of H vary from to 1, and a value of 0.5 implies uncorrelated (random) behavior. As the fBm 
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model indicates, large H values tend to be associated with relatively smooth functions and small H values tend to favor 
highly jagged curves. This feature suggests that H may be useful in quantifying the variability observed in GRB prompt- 
emission light curves. Plotted in Fig. [6] are the extracted //-exponents as histograms for both long and short GRBs. The 
histograms clearly show a displacement in H for the distributions of long and short GRBs, with the short GRBs indicating a 
preference for small values of H. The mean displacement in H raises the interesting possibility of using this feature as a way 
of distinguishing between short and long GRBs. This would be in addition to the currently employed criteria based on T90 
and spectral hardness ratios. Interestingly, the histograms also show a significant overlap in the region of small H exponents 
possibly signaling similarities between the two types of bursts in this range. 

It could be argued that the sizable overlap of the distributions is essentially a consequence of the large dispersion (in 
H) exhibited by both short and long GRB distributions. While it is not known precisely what processes lead to this large 
dispersion in H, we note that the dispersion for the short GRBs is somewhat smaller than the corresponding one for long 
GRBs. If the dispersion is associated with the energetics of the progenitors of the respective systems, i.e., a merger of compact 
objects in the case of short GRBs and the collapse of a rapidly rotating massive star for long GRBs, then one might indeed 
expect a larger dispersion in the //-distribution of long GRBs compared to the corresponding one for short GRBs based 
purely on the difference in the mass range for the respective progenitors. Moreover, additional factors such as the formation 
of an accretion disk, the size of the disk, the mass of the disk, the strength of the magnetic field and the magnitude of the 
accretion rate during the prompt phase, remain largely uncertain. With the added intrinsic variability of the central engine 
itself, we should not be surprised to observe a systematic difference in the extracted Hurst exponents for long and short bursts. 
For completeness, we mention that while the dispersion in H is large for both distributions, the extracted //-value for each 
individual GRB is known reasonably precisely (see Table [I}. 

Another way to examine the //-distrib ution s is to recast the data aga inst the so-called minimum-time-scale parameter, 
MTS, extracted by[ 



MacLachlan ct al 



( 20121 ) and 



MacLachlan ct al 



1 2012al ). Using a method based on wavelets, these authors 



explored the scaling characteristics of GRBs and determined the minimum time scale at which scaling processes dominate 
over random noise processes. Furthermore, the authors have recently shown a direct connection between the extracted MTS 
and the smallest pulse structures extracted by pulse-fitting techniques. In addition to this link with pulses, MTS provides an 
alternate scale (to T90) by which long and short GRBs can be separated. Shown in Fig. [7] are the extracted //-exponents for 
both long and short GRBs versus the MTS (in the observer frame). Short GRBs tend to cluster around small MTS values and 
follow a steep trajectory in the //-MTS plane whereas the long GRBs are distributed over a larger range in MTS and seem to 
follow a gradual power-law- like trajectory. The behavior is a little more clear in panel (b) of Fig. [7| where the MTS is plotted 
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on a log scale: Here the the short and long GRBs indicate a small (~ 30%) positive correlation respectively; the combined 
sample on the other hand shows a larger positive correlation (~ 50%) and an obvious separation of the two distributions with 
MTS. 

Other astrophysical systems for which the Hurst exponent has been extracted includes CVs. These systems, comprising 
tightly-bound binaries (with periods of the order of few hours) and a primary consisting of a compact object (typically a white 
dwarf) and an accretion disk that can accommodate significant mass transfer from the secondary may provide a benchmark for 
gauging the systematics of the extracted Hurst exponents. Indeed, large dispersions in H are found for both optical and X-ray 
light curves of CVs. Interestingly though, CVs apparently tend to favor large ff-exponents i.e., greater than 0.5. This implies 
that the CV distributions are persistent as opposed to a tendency toward antipersistence for GRBs. By their very nature, 
CVs are systems that have built-in periodicity that is readily reflected, in most cases, in the observed emissions from these 
systems. GRBs, on the other hand, are transient phenomena which show very little evidence for periodicities. It's possible 
this simple difference may le ad, in part at least, to the degree of persis t ence or antipersiste nce exhibited by these systems. 



Furthermore, some authors, (jFritz fc Bruch 



1998 



Tamburini et al 



2009 



Anzolin et al 



20101 1. have noted that the extracted 



H's indicate a sensitivity to the strength of the magnetic field of the systems under study, and in particular, the optical and 
x-ray emissions from CVs exhibit different H distributions. Since the optical and x-ray emissions in CVs arise from spatially 
separated regions (the optical from an extended disk and the x-ray in the boundary layer between the inner regions of the 
disk and the surface of the compact object or the polar regions in the case of a highly magnetic system), it is tempting to 
surmise that such a comparison might be fruitful in elucidating the spatial characteristics of GRB jets: Examples include 
the radii and or regions that are commonly associated with the emission sites for prompt gamma-rays (e.g., the photospheric 
radius in the case of a thermal component) and the steeply declining phase of the x-ray light curves (linked with high-latitude 
emission resulting from internal shell collisions). While it is understood that GRBs and CVs are very different systems and 
therefore the translation of the Hurst exponent from one system to the other is likely to be speculative at best, it is intriguing 
nonetheless that a simple scaling parameter may enable us to connect common underlying properties and processes that 
ultimately produce the observed emission in these diverse systems. 



5 CONCLUSIONS 

We have studied the temporal properties of a sample of prompt-emission light curves for short and long-duration GRBs 
detected by the Fermi/GBM mission. By using a technique based on wavelets we have extracted the Hurst exponents for 
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these bursts. This exponent is essentially a scaling parameter that provides a measure of long-range behavior in a time 
series. The physical limits of this index are and 1, where the mid-point (H = 0.5), is an indicator of completely uncorrelated 
(random) processes that contribute to the observed time series. Often times, the H is also associated with the fractal dimension 
(D) of structures by 

D = 2 - H, (14) 

and can be thought of as a measure of the degree of jaggedness of the structures under study. In this sense the H may also be 
indirectly linked to the variability seen in the prompt-emission of many GRBs. Our main results are summarized as follows: 

a) The means of the H distributions for the GRBs in our sample show an offset between short and long GRBs, with the 
short GRBs indicating a preference for smaller Hurst exponents compared to the long GRBs. This offset is potentially an 
independent criterion for distinguishing between long and short-duration bursts. 

b) Compared to short GRBs, long-duration bursts exhibit a larger dispersion in H. The origin of this dispersion is not 
known although it is possible that it is related to the underlying energetics of the different progenitors that produce long and 
short-duration bursts. 

c) No distinct group or clustering is found for H values corresponding to 0.5. This implies that random (or uncorrelated) 
processes, if present, play a lesser role in the production of the observed prompt emission. Moreover, the means of the H- 
distributions for both long and short GRBs indicate a skewness toward values less than 0.5. Overall, this implies that the 
prompt-emission time series exhibit antipersistence. 

Finally, we note that because of the large dispersion in H, there exists a significant region over which the long and short 
bursts overlap. This overlap region raises the interesting possibility of exploring bursts that may possess many more common 
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features than would otherwise be suspected. The case for an intermediate class of GRBs (|Horvathlll998l ; iGao et al.l l2010n 
remains unsettled and warrants further investigation. 
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Table 2. Summary of Long GRBs. 
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Table 3. Summary of Short GRBs. 



GRB Number H SH T 90 [sec] ST 90 [sec] 
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